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We present a general method to study the non-equihbrium behavior of Casimir type fluctuation 
induced forces for classical free scalar field theories, fn particular we analyze the temporal evolution 
of the force towards its equilibrium value when the field dynamics is given by a general class of over 
damped stochastic dynamics (including the model A and model B class). The steady state force 
is also analyzed for systems which have non-equilibrium steady states, for instance where they are 
driven by colored noise. The key to the method is that out of equilibrium force is computed by 
specifying an energy of interaction between the field and the surfaces in the problem. In general we 
find that there is a mapping of the dynamical problem onto a corresponding static one, and in the 
case where the latter can be solved the full dynamical behavior of the force can be extracted. The 
method is used how to compute the non-equilibrium Casimir force induced between two parallel 
plates by a fluctuating field, in the cases of Dirichlet, Neumann and mixed boundary conditions. 
Various other examples, such as the fluctuation induced force between inclusions in fluctuating 
media are discussed. 

PACS numbers: 05.70.Ln, 64.60.Ht 



The Casimir effect arises when the fluctuations of a quantum field is modified by the presence of surfaces or objects 
placed in the field [TH5j . A so called pseudo or critical Casimir force can also arise for classical fields in the presence of 
thermal fluctuations, and although the physical origin of this pseudo-Casimir force is quite different to the quantum 
Casimir effect the two forces have a similar origin from a mathematical stand point. The terminology used to describe 
these types of Casimir forces varies considerably in the literature. We will be interested in the classical, non-quantum 
(or zero Matsubara frequency) Casimir interaction which is dominant at high temperatures and large separations 
between interacting surfaces. The interaction arises in the context of a statistical fleld theory rather than a quantum 
fleld theory and is often referred to as the thermal Casimir effect as in this regime the force is proportional to the 
temperature. When the statistical fleld theory describes a critical point, then one can talk precisely of a critical 
Casimir effect. However even away from a critical point a Casimir type interaction occurs, but this interaction is 
screened rather being long-range. The interaction induced between surfaces or objects can be considered to be due 
to the imposition of boundary conditions on the fleld or due to an energy of interaction with the field. The simplest 
system where the thermal Casimir (or pseudo-Casimir) effect arises is the free scalar field theory, and long range 
Casimir forces arise when the field theory is massless, i.e. where the Hamiltonian of the system is given by: 



The above field theory is purely classical and is, for instance, a simple model for the elastic energy of a surface if the 
fleld 4> represents the height of the surface. If one specifles the boundary conditions for the field on two plates {e.g. 
Dirichlet, Neumann or Robin) the equilibrium thermal Casimir interaction between the plates can be computed from 
the free energy The fluctuating field <j) also describes the order parameter for critical systems, such as binary 

liquids at the critical point [B] , and the occurrence of such a critical Casimir interaction has recently been confirmed 
experimentally [7]. The effect of the Casimir force on the wetting and thinning of '^He and '^He He mixtures has 
also been extensively studied, theoretically [HI H] , numerically [TUl [TT] and experimentally [T2J [13] . A recent review 
of results on the critical Casimir force can be found in [T^. Also the field 4> can represent the phase of the complex 
order parameter in a superfluid state, such as that occurring for '^He [15 , giving an additional contribution to the 
Casmir-like forces in wetting films of "^He. We note that in general critical systems are described by interacting 
field theories, the approach in this paper only applies to free field theories and can only be applied to interacting 
field theories as a Gaussian approximation. Free vectorial field theories describing liquid crystal systems also exhibit 
Casimir type interactions between the surfaces confining the system |16j . In what follows we will study classical fleld 
theories where the dynamics of the flelds are driven by thermal fluctuations or other stochastic noise, i.e. we will 
consider a range of dissipative dynamics for the system. For brevity we shall refer to the forces as Casimir forces 
rather than pseudo Casimir or critical Casimir forces, and at no point will we consider quantum effects. 



I. INTRODUCTION 
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When classical fields are out of thermal equilibrium, one expects significant changes in the Casimir forces from 
those computed at equilibrium. One of the principal problems when analyzing the out of equilibrium Casimir effect, 
say for parallel plates, is to obtain an expression for the force between the two plates which is valid out of equilibrium. 
In previous studies, the stress tensor has been used to study both the dynamical behavior of the force [TTHTn] and the 
force fluctuations in equilibrium j21j. However, only the force in equilibrium can be strictly computed using the stress 
tensor (see the later discussion), therefore there is no general proof that the stress tensor can be used to compute 
forces out of equilibrium for general dissipative dynamics. Results using the stress tensor may, however, be reliable 
for situations close to equilibrium [191 120j . For the quantum Casimir interaction the field dynamics is Lagrangian 
and progress is being made on its out of equilibrium behavior [S]. Another approach to compute out of equilibrium 
Casimir forces is to construct a model with a specified non equilibrium dynamics and to specify by hand the force 
at the wall. For example, in [531 HSIj the dynamical field was related to a particle density and the local pressure on 
the wall is then given by the ideal gas form via kinetic reasoning. We note that in a number of papers where the 
stress tensor is used out-off equilibrium or where kinetic arguments are used to compute the force in parallel plate 
geometries it has been found that the forces exerted on the plates are not equal and opposite. In some of these works 
^} this fact is interpreted as a violation of the Newton action reaction principle. Given these intriguing and 
interesting results it is therefore useful to develop an analysis where we can write down the instantaneous force on 
the plates or surfaces in terms of the unaveraged field variables without making any further assumptions to evaluate 
the force.. 

To achieve this goal, in this paper rather than imposing boundary conditions on the field we specify its energy of 
interaction with the surfaces in the system (for the static problem see |M| for example). In this way we can write 
down the instantaneous force on the wall unambiguously as all forces in the problem are generated by a potential. 
Recently a variant of this approach was also applied to compute dynamical drag forces in fiuctuating classical fields 
pS] This paper presents the full details of the calculations presented in an earlier letter [25] and in addition clarifies 
how our formalism for computing the force is related to the stress tensor. We also discuss the out of equilibrium force 
between plates with mixed boundary conditions for parallel plane geometries. By using the pairwise approximation 
we also show how this method can be applied to compute the out of equilibrium fluctuation induced force between 
small inclusions in fluctuating media. 



II. GENERAL FORMALISM 



A. Energetic formulation of the boundary interaction 



We commence by considering the most general case of a free field theory where the Hamiltonian can be written in 
terms of a general quadratic Hamiltonian |27j . 

H^IJ dxrfx>(x)A(x, x', /)(/.(x'), (2) 

where A is a self-adjoint operator i.e. A(x, x') — A(x',x). Here I represents any suitable free parameter in the 
problem but for concreteness it could be the position of a plate which interacts with the field. For instance one could 
chose we choose 



A(x, x', /) = - [V • k(z, OV - (5(z)ci - S{z - l)c2] 5(x - x'), (3) 

where when ci and C2 are positive this corresponds to a free field theory where the fluctuations of the field (f) are 
suppressed on two plates, one at z = and the other at z = I. The term k(z, /) is a spatially varying elastic constant 
for the field, for instance one could have one value within the two plates and another outside. The induced boundary 
conditions for this theory at each plate are of the Robin form: 

n{0+,l&o+-n{0-,l&o-=cim (4) 
oz oz 

at z = and 

=cMl), (5) 

oz oz 

at z = Z and where the superscripts indicate being infinitesimally to the right and left of the point x. Clearly, 
when K is constant and in the limit where c — > oo one will obtain Dirichlet boundary conditions on the two plates. 
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The instantaneous generalized force acting on the plate at z = Z is thus given by 

Fi = -'^^~-j dxdx>(x)|^A(x,x',O0(x'). (6) 

This is the strict definition as defined by the principle of virtual work, it is valid for any configuration of the field 
4> and position of the interacting surface or object. Clearly simply that the potential energy of a physical system is 
described by a potential V^(x) then the instantaneous force in the direction i for any configuration of the coordinates 

- 

Ji o ' 

OXi 

we are thus simply applying this result to functional potential energies. The equilibrium value of this force is can also 
be written in the familiar form 

{F,)=T^^HZiA)) (7) 

where 

Z{A)= [ d[</)]exp(-/3i7[A]) (8) 



with H as defined in Eq. ([2| and where /3 = 1/T is the inverse temperature. Alternatively the force can be expressed 
using Eq. ^ to obtain 



|a(x,x',0 



A-i(x,x',0, (9) 



where we have simply used the fact that the static correlation function is given by 

(0(x)0(x')>=rA-^(x,x'). (10) 

Before using the above idea to compute the Casimir force out of equilibrium we will briefly show how the energetic 
formulation can be used to derive the standard expression for the average value of the force at thermal equilibrium. 
This is a useful exercise as it shows that our method recovers the stress tensor result in equilibrium and it is also a 
useful reminder as to when the use of the stress tensor is valid. The more traditional derivation in the context of 
quantum field theory can be found for example in |28| . 

Consider the Hamiltonian of a scalar field theory H with Hamiltonian density H such that 

H = J d^iH (11) 

and where the Hamiltonian density consists of a bulk part denoted by Hq and an interaction term with a surface S 
denoted by Hs: 

n = no + ns. (12) 

Here we assume that Hq has a standard quadratic kinetic term [V0]^ plus an interacting term depending only on ip. 
For instance if the surface is perpendicular to the direction i and is at Xi ~ li then we can write 

ns = 6ix,~k)Viq}). (13) 

The force on the surface in the direction i is then given by 

Fi = -^rr- = - / dx — = - / rfx— , (14) 



dk J dli Jvs dk ' 
where the volume Vs in the last integral represents any volume (which can be infinitesimally small) containing the 



surface, this is because of the localized nature of the surface interaction given in Eq. (13). 

In order to make the connection between our energetic formalism and the stress tensor we compute the derivative 
of the Hamiltonian density H in the direction i which can be written as 



(15) 
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we can now use this is Eq. ( 14 ) to give 



Vs 



dn 



Now integrating the last term by parts we find that 

/ dn 



This can now be written as 
where 



Vs 



dn „ f dn 



Vs 



Vs 



dn 



- V, 



dn 

dV,d 



dn 

dVd 



(16) 



(17) 



(18) 



(19) 



is the standard stress tensor. We see that in general there is an additional term in the general force given by the 
second integral above but which clearly must be zero in equilibrium. To see how this vanishes at equilibrium we note 
that the argument of the second integral above is in fact the functional derivative of the Hamiltonian H, i.e. 



dn 



dn 

dVid 



6H 



(20) 



Note that in equilibrium the classical field obeys SH/d(f) — 0, however strictly speaking to show the vanishing of the 
second integral in equilibrium we use the Schwinger-Dyson equation 



(x) 



(0(y)exp(-/3i/))=O, 



as the functional integral is an exact derivative. This can now be rearranged to give 



(0(y) 



SH 



1 



<5(x-y) 



(21) 



(22) 



'(50(x)' /3 

where the angle brackets indicate averaging with respect to the Gibbs-Boltzmann weight exp{—/3H). This then gives 



SU 1 



(23) 



where the last step can be justified by thinking of the Dirac delta function as the limit of a Gaussian. We thus recover 
the equilibrium result 



{F^) 



Vs 



ffx VjTij) 



(24) 



B. Dynamics 

We now consider the dynamical problem where the system in prepared in a state (f> — at the time t — (this could 
have been by cooling the system to a very low temperature for instance) and then letting it relax at some non-zero 
temperature T. We will consider case of a general relaxational dynamics, where the evolution of the field is given by 



dt 



- - J dx' RA{y:,x')(j){x')+r]{x,t) (25) 



where RA indicates the composed operator i?A(x, x') = J (iyi?(x, y)A(y, x'). To satisfy detailed balance with noise 
that is uncorrelated in time we chose the noise correlation to be 



(r/(x, t)f]{x', t')) = 2TS{t - t')R{x, x'). 



(26) 



For the case where i?(x, x') — (5(x — x') we recover the case of nonconserved model A dynamics and when i?(x, x') = 
— V^5(x — x') we have the case of conserved model B dynamics. In what follows the calculation is valid for any self- 
adjoint operator R. 

The formal solution to this equation, for flat initial configuration of the field, is 



(j){x,t)^ / dsdyexp(-(t- s)i?A) (x,y)?7(y,s) 
Jo 

This means that the equal time correlation function of the field is given by, in explicit non-operator notation, 
(0(x,<)0(x',t)) =C(x,x',i) 

= 2T / dsdydy'exp(-(i-s)i?A)(x,y)exp(-(i-s)i?A)(x',y')i?(y,y') 
Jo 

Now we use the fact that the operators A and R are self adjoint to write 

exp {-{t - s)RA) (x'.y') = exp {-{t - s)AR) (y', x') 
This thus enables us to write in operator notation that 

C{t)^2T [ dsexp(-(t-s)i?A)i?exp(-(t-s)Ai?) 
Jo 

Now if we expand the exponential operators in the integral we find 

exp(-(t-s)i?A)i?exp(-(i-s)Ai?) = V -^--(-(t - s))"'+"(i?A)"i?(Ai?)"' 

^-^ n'.m] 

n.m 

^ y-^—(-(t-s))™+"i?(AE)"(Ai?)™ 
^-^ n\m\ 

= Rexp{-{t~ s)2AR) 

The time integration can now be carried out to yield 

C{t) = TA-^ [1 - exp{-2tAR)] 
Now if we Laplace transform this equation (defining Cf[s) — dt exp(— s<)/(t))we find that 



(27) 



(28) 



(29) 



(30) 



CC{s) = 



A- 



(31) 
(32) 

(33) 



Now using this and Eq. (|6| we find that the Laplace transform for the average value generalized force is given by 



mis)) 



T 



dxdx' 



A(x,x',Z) 



sR- 



(x,x') 



However as R does not depend on I we may write 



(£i^,(s)) = / dxdx' 



— (A- 

dv 



sR- 



-1 1 



(x,x') 



sR- 



n -1 



(x,x') 



Now using the equivalence between using Eq, ([7| and Eq. (|9| we may write 



mis))^-§-Hz{A,)), 

s ol 



(34) 



(35) 



(36) 
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FIG. 1: For Dirichlet-Dirichlet boundary conditions the approach to equilibrium for various values of the plate separation 
I (top to bottom I = 1,2,4) for d = 3. Nonequilibrium Casimir force F, obtained by the direct numerical inverse Laplace 
transform of Eq. (38 1, plotted in units of the equilibrium Casimir force {Feq) as a function of time. 

where the operator As is given by 



A, = A + -i?-i 



(37) 



This result is quite remarkable - it means that the Laplace transform of the time dependent Casimir force considered 
here is given by a static Casimir force for another free field theory. It is clear that the result is also valid for the 
force on any surface in the system which interacts with the field. Providing the static partition function is known for 
the corresponding static problem, the corresponding time dependent force can be extracted by inverting the Laplace 
transform. If one takes the limit s ^ in Eq. (361, we find that the static result is recovered from the pole at s = 



III. MODEL A DYNAMICS 



In this section we will analyze the case of model A dynamics i.e. where the dynamical operator i?(x— x') = (5(x— x'). 
This is the easiest case to analyze and it is the case that has been most studied in the literature via the other approaches 
mentioned in the Introduction. 



A. Parallel Plate Geometries 



We now first turn to the case where the imposed boundary conditions are Dirichlet (DD) and the two plates are 
immersed in the fluctuating medium (hence there is fluctuating medium on both sides of each plate). We take the 
total length of the system to be L which is fixed, and place the plates of area A a distance I apart. Standard results 
on the screened Casimir interaction [3l lU [3^ give 

9 AT 

mis)) = - _ _ /(^, .), (38) 

.s(47r) — r(V) 

where d is the dimension of the space, A is the area of the plates and. 



1 - exp(-2ly'fc2 + f ) 
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The equilibrium behavior is easily extracted by examining the pole at s = which yields, as anticipated, the standard 
equilibrium Casimir force 



ATr{d)Cid) 



where T is Euler's gamma function and ^ is the Riemann zeta function |31) 



Cid) 



OO 

V — 

n=l 



(40) 



(41) 



The full time dependence of the force, starting at zero at t = and relaxing to the equilibrium value above, can be 
extracted by direct Laplace inversion of Eq. ( 38 1 . Fig. ([T ) shows the approach to equilibrium for three different plate 
separations. Clearly, the relaxation times increase with plate separation and this is due to the fact that the underlying 
dynamics is diffusive and hence P sets a time scale. This is also clearly evident in Fig. ^ that shows the collapsed 
rescaled force curves obtained by plotting with the time units rescaled by P. 

Useful analytic expressions for the early and late time behavior of the non-equilibrium force can also be obtained 
from Eq. ( 38 1 , because of the technical nature of their derivation they are relegated to the Appendix. We find that 



the temporal derivative of the out of equilibrium force is given by 



' dt ' 



2AT d 



OO 



The above expression may also be written in the form 



cxp(-^). 



d_F^ 
^ dt ' 



AT 



AT 



d 



2(87r)ir 



2(8^)- 



r^t^ dl 



Z^cxp( j^) 



(42) 



(43) 



Clearly because the underlying dynamics is diffusive P sets a time scale. The short term behavior of the force, 
t/P ^ 1, can be obtained directly from Eq. (42) and is given by 



The long time asymptotics, t/P ^ 1, follow directly from Eq. (43) and are given by 



(Fi) ^ {Fi),, + 



AT 



d{8Tr)Hi 



(44) 



(45) 



An interesting thing about Eq. ( 43 ) is that one sees explicitly the appearance of the eigenvalues for Dirichlet boundary 



conditions in the sum on the right hand side. The first term can be seen to be due to the bulk on the exterior of the 
system. This can shown be by taking the limit I — > oo and expressing the sum as a Riemann integral gives the left 
hand side to be equal to zero, i.e as one would expect there is no force. Of course for a system where I — !■ oo the 
system is always out of equilibrium. 

The agreement between our asymptotic expressions and exact results obtained by numerical inversion of Eq. ( 38 ) 
are shown in Fig. ([2]). As mentioned above the late time correction is independent of I. This is because the medium 
between the two plates has a relaxation time t{1) ~ P whereas the slowest relaxation times in the system 
are associated with the medium outside the two plates and hence at late times the correction is dominated by the 
relaxation of the external system in the thermodynamic limit L ^ oo. This diffusive relaxation is responsible for the 
power law approach to equilibrium. The above result is also valid for Neumann-Neumann (NN) boundary conditions 
as the static screened problems have the same force for both DD and NN boundary conditions. We note that this 
problem has been studied in |20| assuming that the stress tensor can be used to compute the force out of equilibrium. 
Our results for DD boundary conditions agree with that given by the stress tensor, though no explicit formula for 
the force is given in [5D] we are able to show that this should give the same result for the force (see later). Indeed 
the numerical curve given for the DD case in [20] closely resembles ours in Fig. (jlj and the asymptotic form of the 
late time decay to the equilibrium force agrees with ours. It is also clear that within this energetic formulation that 
the force on the two planes, for any boundary conditions, is equal and opposite as it is given via the equilibrium force 
of another free field theory. In [20 the difference between the DD and NN boundary conditions can be seen to be 
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FIG. 2: Nonequilibrium Casimir force, F, in units of the equilibrium Casimir force {Feq) at late times and early times (inset) 
for d = 3. Symbols were obtained by numerical inverse Laplace transform of Eq. (38|. Each plot has only one set of symbols 



because curves for different I values collapse. Dashed lines are the corresponding approximations from (Eqs. ( |44[ ) and (|45[)). 

due to the presence, for NN boundary conditions, of the zero (constant) mode in the z direction (perpendicular to 
the plates) in the computation using the stress tensor. This is the only difference between the two results, however 
it is not clear to us how the zero mode can influence the force as it does not see {i.e. it is not effected by) the NN 
boundary conditions. 

To see whether, in general, a stress tensor computation agrees with our result, we can exploit the fact that the right 
hand side of Eq. ( 36 ) is for an equilibrium force, for the Hamiltonian 



(46) 



the subscript st is to make clear that the static or equilibrium measure with Hamiltonian ( 46 ) is used to compute 



observables. Now away from (but still possibly infinitesimally close to) the boundary the energy density is given by 



The stress tensor for this theory is given by 



V,(/)Vj0, 



(47) 

(48) 
(49) 



and the average force on any volume V is given by the average of integral over the bounding surface S of this stress 
tensor 

(50) 



Js 



The Laplace transform of the force on the plates is thus given from Eq. ( 36 ) by 

s Js 



^3 I St 



(51) 
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where, to avoid possible confusion the subscript dy indicates averaging over the noise for the dynamical problem and 
S indicates the surface of the plates. However we can now use the relation Eq. (33) to write 

£((^(x,O0(x',i))d,)(.s) = i(0(x)0(x')).t, (52) 

and we may thus write 

{r.A^))st = {cm<i>f),y)is) + '-cm,y)is)) /:((v.0v,0),,)(.). (53) 

Now we use the fact that if 0(x, 0) = then {(f>{x, 0)(j){x' , 0))dy = and thus 

'-Cm,y){s) = l£((^02^,^)(.) (54) 

The inverse Laplace transform of the static stress tensor 7ij{s) can thus be treated as an effective dynamical stress 
tensor 

{Tffit)),y = C-'(%,{s)),,,it) (55) 

with 

T^tW = % ( m' + —]- V.0V,</' (56) 



2 dt _ 

and (f> the dynamical field. Therefore computing forces with this dynamical stress tensor will give the same forces 
as those given via our boundary energy derivation. Let us emphasize here that the effective dynamical stress tensor 
written here is by no means a universal one, it depends on the precise dynamics of the system and the choice of initial 
conditions. 

We see immediately from Eq. ( 56 1 that at late time the average value of the time derivative term will go to zero 
and we will recover the standard form of the stress tensor . We also see that when computing the force on a plate 
with Dirichlet boundary conditions the temporal derivative does not contribute as </> is zero at the surface. This is 
why the results of [2D] agree with ours for DD boundary conditions. 

In the case of NN boundary conditions it can be shown that the temporal derivative term cancels out the contribution 
from the zero mode and thus gives exactly the same force as for DD boundary conditions. Similarly for DN boundary 
conditions this temporal derivative term ensures that the force at the D boundary is equal in magnitude but opposite 
to that at the N boundary, this is in contradiction with the corresponding result given in [50] where the forces at the 
boundaries are not equal and opposite. 

In the case of DN boundary conditions the screened static result is 

9 AT 
s(47r)T-r(^) 

but here, 

,d-2,, ypTfexp(-2/ypTi) 



I(l,s)^- / k^-^dk^ ^ — — / " , (58) 

7 l + exp(-2;yfc^) ^ ' 

The equilibrium behavior is easily extracted by examining the pole at s = which yields, as anticipated, the standard 
equilibrium Casimir force for DN boundary conditions 

^^'^^^-(i6.)^r(^)/'^' ^''^ 



where 



C* W = E = (1 - 2i-'^)C(d) (60) 



n=l 
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is positive for d > 1 and so the force in this case is repulsive. Exactly the same analysis as that given above can be 
applied to give 



,dF, 



2^2 



Pn 



I 2AT d ^ 1 ,x„-i / « ^ 



From this we obtain immediately the short time asymptotic behavior of the force 

(Flit)) 



2AT , P , 



(61) 



(62) 



we see it is thus opposite that of the case of DD and NN boundary conditions. The large time asymptotic behavior 
of the force can also be extracted by using another form of the Poisson summation formula, namely 



l + 2^(-l)"exp(- 



2t 



-) = 



/2'Kt 

IT 



1^ exp( 



which in the late time limit gives 



and thus 



oo 

E 



' dt ' 



AT 

2(8^)5r 



and hence 



{Fi) ^ {Fi),, + 



AT 



(63) 



(64) 



(65) 



(66) 



Therefore we see that the time correction to the equilibrium result for DN boundary conditions is exactly the same 
as that for the DD and NN case. However this means that the intermediate force must overshoot its equilibrium 
value. Fig ([s]) shows the full time dependent approach to equilibrium for three different plate separations obtained by 
direct numerical inversion of Eq.(57). The overshoots at intermediate times are clearly visible. This agrees with the 



result of [20] for the force calculated at the Dirichlet wall as it should from our arguments above stating that one can, 
as was done in [20|, use the usual expression for the stress tensor at the plate with Dirichlet boundary conditions. 
This overshoot may possibly be explained in the following way, the effect of the external medium for both types of 
boundary conditions is to cause an additional, temporally decaying, repulsion between the plates, this is irrespective 
of whether the equilibrium force between the plates is repulsive or attractive. Thus in the case where the equilibrium 
force is repulsive we will clearly have an overshoot effect. This picture is backed up by the fact that the decay is 
independent of the distance between the plates. 



B. Effect of Temperature Changes 

The method can also be used to examine the dynamics resulting from a sudden change in temperature, from say 
Tq where the system is in equilibrium to a temperature T. In this case, the initial configuration of the field (/'(x, 0) 
has the correlation function 

((/.(x,O)0(x',O)) =roA-i(x,x',0. (67) 



Solving the equation of motion Eq. ( 25 1 with this initial condition yields the time dependent correlation function 

C(x, x', t) = To exp(-2tA)A"i + T{1 - exp(-2tA))A~i (68) 

and where the second term is exactly the same as that arising for flat initial conditions. Taking the Laplace transform 
of this gives s£C(x,x',i)(s) = Tq [A-i(x,x,r) - Aj^{x,x,l')] +TAj^{x,yi,r). and we find 

mis)) ^ ^|ln(Z(A)) + ^-^|ln(Z(A,)) (69) 
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FIG. 3: For Dirichlet-Neumann boundary conditions, the approach to equilibrium for various values of the plate separation / 
for d = 3. Nonequilibrium Casimir force F, obtained by the direct numerical inverse Laplace transform of Eq. ( |57| l, plotted in 
units of the equilibrium Casimir force [Feq] as a function of time. 

where we have used the fact that ^ ln(Z(As)) is independent of the temperature. For DD boundary conditions this 
gives the Umiting behaviors 

{Fl{t)) ^ (fL)e,To- ^^5.T-^°^ "^P(-|) T»l 

d[87r)2t2 t 



(70) 



IV. NON-EQUILIBRIUM STEADY STATE FORCE - MODEL A FORCING WITH COLORED NOISE 



One can also consider the behavior of the Casimir force for relaxational dynamics where the deterministic 
forcing term is of model A type but where the forcing noise is colored in time such that (r/(x, t)?7(x', t)) ~ 
TS{x-x')u}exp{-uj\t-t\') i.e 



dt 



'1^ 



V 



(71) 



Here T represents an energy scale, a; a frequency and the resulting steady state is not an equilibrium one. The 
average value of the force in the steady state regime can be computed using the same formalism above and we find 
that the correlation function of the field is given by 



which yields 



C(x, x', 0.) = T [A(x, x', - A2„(x, x', 0"'] 



{Fiiio)) = -\j dxdx'|A(x,x',OC(x,x',L.) 
= r|[ln(Z(A))-ln(Z(A2.)] 



(72) 



(73) 



Hence again we find that one can compute a force in a non equilibrium system from knowledge of static screened 
systems [32] . Note that in the limit — )■ oo we recover the white noise equilibrium result of Eq. ([t]) . Fig. Q shows the 
frequency dependence of the non-equilibrium force obtained from Eq.(73l for a two plate system with DD boundary 
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FIG. 4: Steady state pseudo-Casimir force, F, for colored noise with Dirichlet-Dirichlet boundary conditions in units of the 
equilibrium Casimir force {F^q) as a function of ijj (in units of for I = 1 and d = 3. Inset: Force per unit area, F (in units 
of kBT/lf), where Iq is the unit of length), as a function of plate separation I (in units of Iq) for uj = 1 ((in units of Iq'^)) and 
d — 3 (solid line). The equilibrium Casimir force (dashed line) is shown for comparison. 



conditions. This result is also that same as that for NN boundary conditions by the equivalence of the corresponding 
static problems. Again P sets a timescale and we see that for lu ^ the force, F tends to the equilibrium white 
noise value, F^q, as expected, while for ui ^ l^^, F <C Feq and a.s lu 0, F vanishes. The inset to Fig. (jZJ) shows how 
the force depends on plate separation for fixed uj. Again, equilibrium behavior is recovered for large I { uj ^ ^^^), while 
for small plate separations the force changes qualitatively scaling as l^^. We note that the result Eq.(73) agrees with 
a computation for the same system where the steady state Casimir force was computed using the stress tensor |17) . 
Fig. Q shows the frequency dependence of the non-equilibrium force obtained from Eq. ( 73 ) for a two plate system 
with DN boundary conditions. One can see that the qualitative behavior is the same as for DD boundary conditions, 
tending to zero and the equilibrium white noise value for small and large lo respectively. In contrast, the scaling of F 
for small uj is different. This is most clearly manifested if we look at how the force depends on the plate separation 
at fixed u) as shown in the inset to Fig. ([5]). While equilibrium behavior is recovered for large I {uj ^ l^'^), for small 
plate separations the force changes qualitatively becoming almost insensitive to the plate separation. This difference 
is again highlighted in Fig.([6| which shows an explicit comparison between the curves in Fig. Q and Fig. ([5|. 



V. PAIRWISE APPROXIMATION FOR SMALL DEFECT REGIONS 



We can also study the force between two small defect regions in an elastic fluctuating medium using the pairwise 
approximation which neglects n-body effects. For instance one could have two small volumes Vi and V2 separated by 
a distance I in which the held acquires a mass. This will give a Hamiltonian of the form 

H^l fd^ [V0(x)]2 + I f dxcf,\^) + ^ f dx 02 (x) (74) 

The pairwise approximation is equivalent to evaluating the partition function for the Hamiltonian above to second 
order in the cumulant expansion. We find that the effective equilibrium potential arising between the two regions, 
when / is much greater than their sizes, is 

nO = -^^^^i^Gg(0, (75) 
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FIG. 5; Steady state pseudo-Casimir force, F, for colored noise with Dirichlet-Neumann boundary conditions in units of the 
equilibrium Casimir force (Feq) as a function of uj (in units of for / = 1 and d = 3. Inset: Force per unit area, F (in 
units of kBT/lg, where Ig is the unit of length), for Dirichlet-Neumann boundary conditions, as a function of plate separation 
I (in units of lo) for lj = 1 ((in units of Ig'^)) and d = 3 (solid line). The equilibrium Casimir force (dashed line) is shown for 
comparison. 




FIG. 6: Explicit comparison of the steady state pseudo-Casimir force, F, for colored noise with Dirichlet-Neumann (red solid 
line) and Dirichlet-Dirichlet (blue dashed line) boundary conditions in units of the equilibrium Casimir force (Feq) as a function 
of u! (in units of /~^) for / = 1 and d = 3 . Inset : Comparison of the absolute value of the force per unit area, F (in units 
of ksT/lQ, where lo is the unit of length), for Dirichlet-Neumann (red solid line) and Dirichlet-Dirichlet (blue dashed line) 
boundary conditions, as a function of plate separation I (in units of Iq) for cj = 1 ((in units of Iq^)) and d = 3. 



where Go = — V~ is the unscreened Coulomb potential in d dimensions. Note that in Eq. (75) it is the function Go 
evaluated at I vifhich is squared, not the operator . Now we consider how this force evolves towards its static value 
from the initial conditions where = throughout the system. Applying the theory developed above we find that, 
in this same pairwise approximation, the time dependent force between the two defects is given as the derivative of a 
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time dependent potential V{l,t) whose Laplace transform is given by 

mi,s) = -^^^^^^Glil), (76) 

where Gs — (— V'^ + f)^^- integral representation can be found for arbitrary dimension d, however the result in 
d — 3 takes the particularly simple form 

^rn +^ TC1C2V1V2 r r i s /^^n 

where erfc is the complementary error function. The effective interaction between the above types of defects for 
colored driving noise can also be derived via Eq. (73 1. Within the pairwise approximation we find that the steady 
state force in this case is obtained from the effective potential 

m-) = -^^^^[G^(0"GL(/)]. (78) 

^^^'""^ ^ - ^32^2^2^' - exp(-2V^0)- (79) 

Hence at large separations Plu ^ 1 the interaction is the same as the that for white noise but when Pco ^ 1 the 
effective potential behaves as —y/uj/l. 



and for d — 3, for example, we find 



VI. CONCLUSION 



In this paper we have studied how the Casimir interaction due to a fluctuating scalar field between two plates 
behaves as a function of time in the case of relaxational dynamics, driven by white noise and obeying detailed 
balance. In particular we have shown how it evolves toward its equilibrium value from an initial state where all the 
field fluctuations are suppressed. We have also analyzed the steady state behavior of the force induced when the field 
dynamics does not obey detailed balance, notably we have analyzed what happens when the deterministic part of the 
dynamics is relaxational but the noise is colored. Previous studies on the dynamical Casimir effect concentrated on 
steady state non equilibrium dynamics or dynamics close to equilibrium and considered Dirichlet boundary conditions 
assuming that the equilibrium stress tensor could be applied to compute the force. Our formalism marks a major 
advance that overcomes these restrictions by allowing the time-dependent force to be evaluated unambiguously via 
an expression for the energy of the field. The method presented is very general and it would be interesting to analyze 
the temporal behavior of the Casimir force for other types of dynamics. We have restricted ourselves to model A 
type dynamics as the resulting static results necessary to extract the temporal behavior of the force are known. The 
interested reader will see that, for instance, in the case of model B dynamics one must know how to compute the 
Casimir force with a non-local Gaussian action. In addition one must also determine whether the model B dynamics 
conserves the order parameter within the plates or conserves it globally. Our method applies to the case where the 
dynamics is globally conserved as the operator R is assumed to be independent of the plate positions. The interested 
reader can find a discussion of these points in '33' . Another extension of the results here would be to actions with 
higher order derivative terms such as the Helfrich action for membrane fluctuations |34| . The main problem here is that 
the corresponding static results are very complicated even in the case of planar geometries |35j . However the method 
should be relatively straightforward to apply in the context of the pairwise approximation for the interaction between 
membrane inclusions j36H38| which can be used to analyze the interaction of defect regions of different elasticity 
and bending rigidity in membranes. In addition the method used here can be applied to Brownian hydrodynamical 
dynamics which is the relevant dynamics for membrane fluctuations [39j . 

This research was supported in part by the National Science Foundation under Grant No.PHY05-51164 (while 
at the KITP UCSB program The theory and practice of fluctuation induced interactions 2008). DSD acknowledges 
support from the Institut Universitaire de France. AG acknowledges support from a James S. McDonnell Foundation 
Award. 



Appendix A: Analysis of time dependent force for parallel plate geometry with model A dynamics 



Here we carry out an analytical inversion of the Laplace transforms in Eq. ( 38 ) for DD boundary conditions parallel 
plate geometry under model A dynamics. We will use a number of standard textbook properties of Laplace transforms 
which can be found for instance in [29\. The computation for ND computations follow with only minor variations. 
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Starting with the DD case, using the fact that F;(0) = we may write Eq. (38) as 



Now we can use the result 



V2(47r)^r(^^) 

£[exp{-pt)f{t)]{s) = Cf{s+p) 



to obtain 



v2(47r) 2 1 1^:^) J 



y2(4^)^r(^) 

where the Laplace transform of fn{i) is given by 

C-fn{s) = -y/sexp((— \/2Zn-v/s). 

We now note that 

£[crfc(-^)](s) = - exp(— 2a\/s) 
where erfc denotes the complementary error function defined as 



Eq. (A5) can be written as 



which thus gives 



which means that 



2 

erfc(2;) = —j= / du exp(— u^) 

Jz 



^da^^m^'^'^Ti^^^'^ ^ yiexp(-2ayi) 



C ^ y/s exp{—2ay/s) — cxp{ 

\/n t 



"2tl +i§ 



1 /2„2 

/„(t) = -^cxp(-— ; 



1 Pn' 



2^5 2ti 



Putting this together and inverting the Laplace transform yields 

, dFi , AT 



dt 



3 Nil,J±+2 



2,n2 



exp( 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(AlO) 



which gives Eq. (42 1. . Short time correspond to the regime where /t^ 1 and hence in this regime the dominant 
behavior is 



' dt ' 



ATP , l\ 

— 5-i^+rexp(- — ). 



B7r)2t 2 



2t' 



(All) 



Thus at very short times we find that 

The long time asymptotics can be obtained by using the Poisson summation formula 



1 + 2 ^ exp(- 



2^2 



/27rt 



2t 



27r2n2/ 



l + 2^exp( —) 



(A12) 



Using this we find that the time derivative of the average force can be written as Eq. ( 43 ) 
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